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Abstract 

Exclusive double diffractive (EDD) Higgs boson production is analyzed in the frame- 
work of the Regge-eikonal approach. Total and differential cross-sections for the 
process p + p—>p + H + p are calculated. Experimental possibilities to find Higgs 
boson at LHC are discussed. 
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1 Introduction 



LHC collaborations aimed at working in low (TOTEM)- and high (CMS)-pr regimes 
related to typical undulatory (diffractive) and corpuscular (point-like) behaviours of the 
corresponding cross-sections may offer a very exciting possibility to observe an interplay of 
both regimes [1]. In theory the "hard part" can be (hopefully) treated with perturbative 
methods whilst the "soft" one is definitely nonperturbative. 

Below we give an example of such an interplay: exclusive Higgs boson production by 
diffractively scattered protons, i.e. the process p + p^p + H + p, where + means also a 
rapidity gap. 

This process is related to the dominant amplitude of exclusive two-gluon production. 
Driving mechanism of the diffractive processes is the Pomeron. Data on the total cross- 
sections demands unambiguosly for the Pomeron with larger-than-one intercept, thereof 
the need in "unitarisation" . 

As will be seen below the detection of the Higgs boson (in the bb mode) at LHC in 
the double diffractive regime looks fairly possible. 

2 Calculations 

In Figs. 1,2 we illustrate in detail the process p + p^p + H + p. Off-shell proton-gluon 
amplitudes in Fig. 2 are treated by the method developed in Ref. [2], which is based on 
the extension of Regge-eikonal approach, and succesfully used for the description of the 
HERA data [3]. 

The amplitude of the process p + p — > p + H + p consists of two parts (see Fig. 2). The 
"hard" part F is the usual gluon-gluon fusion process calculated in the Standart Model [4]. 
"Soft" amplitudes T 12 are obtained in the Regge-eikonal approach. 

We use the following kinematics, which corresponds to the double Regge limit. It is 
convenient to use light-cone components (+, — ; _L). The components of momenta of the 
hadrons in Fig. 2 are 

Pi 

P2 

Pi 

P2 

q 

Q2 

£1.2 are fractions of protons momenta carried by gluons. For two-dimensional transverse 




(1) 



= (q+,q-, q) , 

= q + Pi - Pi = q + Ai , 

= -g + P2 - p' 2 = -q + A 2 , 
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vectors we use boldface type. From the above notations we can obtain the relations: 
t - A 2 ~ A ^( 1 +6, 2 )+ei> 2 „ 2 

Cl,2 - ^1,2 — \ 7 — -^1,2 > ?1,2 U (2 J 

1 - ?1,2 



COS 0o 



AiA 



2 



|Ai||A 2 



M 2 H - 66s + *i + *2-2Vti*2Cos^ 
(pi + g) 2 ~ m 2 + g 2 + V2sq^ = Si 
(j>2 ~ qf — m 2 + q 2 - V2sq+ = s 2 ■ 

Physical region of diffract ive events with two rapidity gaps is defined by the following 
kinematical cuts: 

0.01 GeV 2 < \t 1>2 \ < 1 GeV 2 , (3) 
M 2 

Un ^ -T^- < 6,2 < Uax = 0.1 , (4) 

(V 3 ^- V^hf < k < (v^+v 7 ^) 2 (5) 
k = 66« - M 2 H <^M 2 H 

We can write the relations in terms of y^ 2 and yu (hadron and Higgs boson rapidities 
correspondingly) : 

6,2 ^ ^e** , (6) 
IVh] < Vo = In 



y ~ 2.5 for = 14 7W , M# = 115 GeV , 

i , i, (i-6, 2 ) 2 ^ n 

^.2 = O 111 2 7 - 9 

2 m 2 — t^ 2 

The calculation of the amplitude is based on the nonfactorized scheme. The contribu- 
tion of the diagram depicted in Fig. 2 is obtained by integrating over all internal loop mo- 
menta. It was shown in [5], that the leading contribution arises from the region of the inte- 
gration, where momentum q is " Glauber-like" , i.e. of the order (k + m 2 / y/s, k_m 2 / y/s, km), 
where k's are of the order 1. The detailed consideration of the loop integral like 

d 4 q f(q,p 1 ,p 2 ,A 1 ,A 2 ) ^ 



(27i) 4 (q 2 + iO)(q 2 1 +iO)(ql + iO) 
shows that the main contribution comes from the poles 

q 2 = V2~s^q. - qi 2 = , 
q 2 = -V2sS 2 q+ - q 2 2 = . 
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In this case 



where 



qi \q), (8) 



q 1 2 =q 2 + A 1 2 + 2|q||A 1 |cos(0+p, 
q 2 2 =q 2 + A 2 2 -2|q||A 2 |cos(0-^) 

In this region \q 2 jM\\ -C 1, and for the ggH-vertex F^ v in the first order of the strong 
coupling we have the usual [4] expression 



r2 a s G F 



F gg ^ H = M H^^m , (io) 



^2 ; + ^2 



1 - i - i0 / I 1 + Jl - - + iO 

v / \ V v 



where 77 = M^/Amf, Gp is the Fermi constant. NLO K-factor 1.5 for the gg — > H process 
is included to the final answer. 

Taking the general form for T-amplitudes that satisfy identities 

q a T° a = 0, tfr% = 0, (11) 

and neglecting terms of the order o(£j), the following expression is found at |^| < 1 GeV 2 : 

/ pn pi \ 

T°(p, q, ft) = [G, a ~ j^ q j T g D p ^ gp ( Si , t t , qq t ) (12) 





= sv 


qi,^qa 
qqt 


pQi 
tl 


= p»- 


pqi 

% 

qqi 


P q -- 
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= Pa~ 


pq 

q%,a 

qqi 



For T®^^ we use the Regge-eikonal approach [2, 6]. At small U it takes the form of the 
Born approximation, i.e. Regge factor: 



2 \ a p(*») 

so - qqi - m 2 



T g UM^,qq>) = c 9P le-^ S ; ™ ™ ) , . (13) 



1 r 2 
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where a P (0) = 1.203, a' P (0) = 0.094 GeV~ 2 , r 2 p = 2.477 GeV" 2 are found in [6]. Pa- 
rameters c gp ~ 3.5, r 2 p = 2.54 GeV -2 are defined from fitting the HERA data on elastic 
J/ty production, which will be published elsewhere. The upper bound for the constant 

~ 2.3(3.3) can be also estimated from the exclusive double diffractive di-jet produc- 
tion at Tevatron, if we take CDF cuts and the upper limit for the exclusive total di-jet 
cross-section [7]. The effective value c gp = 2.3 corresponds to the case, when the Sudakov 
suppression factor is absorbed into the constant, and c gp = 3.3 is obtained when taking 
into account this factor explicitely. 

The full amplitude for Higgs boson production looks as follows: 



, f d*q 8F»"(q 1 ,q 2 )T°(p 1 , q, gi)T£(P2, -q, q 2 ) 

l pP ^ P H P ^J (2?r)4 {q 2 + i0){q 2 + i0){q 2 +i0) ■ ^ 



Factor 8 arises from the colour index contraction. Let I 2 = —q 2 ~ q 2 , y p =< Uh >— 
and contract all the tensor indices, then the integral (14) takes the form 

T pp ^ pHp ~ 4p eb{tl+t2) J^2 [-M*) ■ 8i W " I , (15) 
b = a' P (0) In ( £-) + b , (16) 



L 



M% dl 2 ( / 5 



\ 2a P (0) 



, , ~ 1.88 GeV" 2 , (17) 

l A \s - m 2 + I 2 /2 J ' V ; 

where M H = 100 GeV, and s — m2 — 1 GeV 2 is the scale parameter of the model that is 
used in the global fitting of the data on pp(pp) scattering for on-shell amplitudes [6]. It 
remains fixed in the present calculations. 

If we take into account the emission of virtual " soft" gluons, while prohibiting the real 
ones, that could fill rapidity gaps, it results in the Sudakov-like suppression [8]: 



F s (l 2 ) = exp 



3 M f 4 dp T 2 (M H 2 ~ 
"27 / ^^HW, 



(18) 



and to the new value of the integral (17): 



r M % dl 2 ( I 2 \ 2a p(°) 

'•"L T^ U-m' + Pft ) — 0.38 GeV -2 , M H = 100 GeV. (19, 

In this case the total cross-section becomes 24 times smaller, than without the factor F s . 

Unitarity corrections can be estimated from the elastic pp scattering by the method 
depicted in Fig.l, where 



Tx — Tpp^ P H P , (20) 
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V(s, q r ) = 4s(27r) 2 (5 2 (q T ) +4s J d 2 he 1 ^ 



T% nit -( Pl , p 2 , Ax , A 2 ) 



1 /■ orq T d c( T 



lQss' J (2vr) 2 (27r) 2 
■ , q^) , 

Air = Ai - q T - q T , 
A 2T = A 2 + q T + , 



V(s , q T ) • T x (pi - qr,P2 + qr, A iT , A 2T ) 



and 8pp_p P can be found in Ref. [6]. It reduces the integrated cross-section by the factor 
about 14. 



3 Results and discussions 

We have the following expression for the differential cross-section: 

rlrr -\rpUnit. 12 

aa _ n \ 1 pp^pHp\ 

dt^hd^d^ ~ 8s(27t) 5 v /z A 

A = k 2 + 2(ti +t 2 )K+ (h -t 2 f < 

By partial integrating (21) we obtain the cross-sections da/dt and dojd^. The first 
result of our calculations is depicted in the Fig. 3. The antishrinkage of the diffraction 
peak with increasing Higgs boson mass is the direct consequence of the existence of the 
additional hard scale Mh, which makes the interaction radius smaller. The £ dependence 
is shown in Fig. 4. 

The second result is the total cross-section versus Higgs boson mass. We obtain the 
following results for the total cross-section of the EDD Higgs boson production: 



Cgp 


Mh (GeV) 


Cp+p^p+H+p (fb) 


LHC 


TeVatron 


no Sud. suppr. 


Sud. suppr. 


no Sud. suppr. 


Sud. suppr. 


3.5 


100 -> 500 


110 -> 57 


4.6 -> 0.14 


12 -> 0.4 


0.5 -> 0.001 


2.3(3.3) 


100 -> 500 


20 -> 11 


3.6 -> 0.11 


2.2 -> 0.08 


0.4 -> 0.0009 



(21) 



The above results for LHC energies are depicted in Figs. 5-8. The form of the curve 
originates from the standard gluon-gluon-Higgs vertex and has a peak near M H — 2m t . 
For the case of Sudakov-like suppression the cross-section vanishes faster with Mh- 

It is useful to compare our result with other studies. Results quite close to ours (with 
the normalization to the CDF data, c gp = 3.3) were given in Ref. [8], where the value 
of the total cross-section is about 3 fb. In both cases the most important suppressing 
in the mass region M H > 100 GeV is due to (perturbative) Sudakov factors, while the 
nonperturbative (absorbtive) factors play relatively minor role. 

Results of other authors were considered in details in [9]. The highest cross-section 2 pb 
for Mh = 400 GeV at LHC energies was obtained in Ref. [10]. They used a nonfactorized 
form of the amplitude and a "QCD inspired" model for gp — > gp amplitudes, taking into 
account the nonperturbative proton wave functions. Even if we multiply the result of 
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Ref. [10] by the suppressing factor, it will be larger than ours. This could serve as the 
indication of the role of nonperturbative effects. Our model is based on the Regge-Eikonal 
approach for the amplitudes, which is primordially nonperturbative, normalized to the 
data from HERA on ^p — > J/^p and improved by the CDF data on the exclusive di-jet 
production. 

To estimate the signal to QCD background ratio for bb signal we use the standart 
expression for gg — > bb amplitude and assumptions [11]-[14]: 

• possibility to separate final bb quark jets from gluon jets. If we cannot do it, it will 
increase the background by two orders of magnitude under the 50% efficiency 

• suppression due to the absence of colour-octet bb final states 

• suppression of light fermion pare production, when J z>to t = (see also [15], [16]) 

• cut E T > 50 GeV {6 > 60°), since the cross-section of diffractive bb jet production 
strongly decreases with E T . 

The theoretical result of these numerical estimations is 

Signal (pp — > pHp — > pbbp) 



QCD background 



> 



>9.5.10- 6 |/| 2 ^ Q0 ^, (22) 

where |/| 2 ~ 0.5 -> 3 for Higgs boson masses 100 -> 350 GeV. For M H ~ 115 GeV 

Signal (pp -> pHp -> pbbp) > ^ GeV 

QCD background ~ ' AM ' 1 ' 

where AM is the mass resolution of the detector. Similar result was strictly obtained 
in [11], [12]. More exact estimations of the above ratio by fast Monte-Carlo simulations 
and the total efficiency of EDD Higgs boson production will be publised in a forthcoming 
paper. 



4 Conclusions 

We see from the result that there is a real possibility to detect Higgs boson using the 
usual bb signal under the luminosity greater than 10 32 in EDD events at LHC. Accuracy 
of the mass measurements could be improved by applying the missing mass method [17]. 

The low value of the exclusive Higgs boson production cross-section obtained in this 
paper is mainly due to the Sudakov suppression factor (18), the full validity of which is 
not obvious for us, because the confinement effects can strongly modify the "real gluon 
emission" . 

It is interesting that in spite of different models and quite different ways of account of 
absorbtive effects in our paper and in Ref. [8], the final results (see Fig. 8) appeared to be 
quite close. 
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Certainly, the cross-sections may be several times lager due to still not very well known 
non-perturbative factors. 

It is possible to generalize our approach for the exclusive production of other particles 
like Xco, Xbo, radion, K.-K. gravitons, glueballs. In this case cross-sections can be lager 
than for EDD Higgs boson production, and some other important investigations like 
partial wave analysis and measurements of the diffractive pattern of the interaction could 
be done. 

Aknowledgements 

We are grateful to A. De Roeck, A. Prokudin, A. Rostovtsev, A. Sobol, S. Slabospitsky and 
participants of BLOIS2003 workshop and several CMS meetings for helpful discussions. 
This work is supported by the Russian Foundation for Basic Research, grant no. 02-02- 
16355 

References 

[I] V.A. Petrov, Proc. of the 2nd Int. Symp. "LHC: Physics and Detectors". (Eds. 
A.N. Sissakian and Y.A. Kultchitsky), June 2000, Dubna. P.223; V.A. Petrov, A.V. 
Prokudin, S.M. Troshin, N.E. Tyurin, J. Phys. G: Nucl. Part. Phys. 27 (2001) 2225. 

[2] V.A. Petrov, Proceedings of the Vllth Blois Wokshop (Ed. M. Haguenauer et al, 
Edintions Frontieres; Paris 1995). 

[3] V.A. Petrov and A.V. Prokudin, Phys. Atom. Nucl. 62 (1999) 1562. 

[4] A. Kniehl, Phys. Rept. 240 (1994) 211. 

[5] A. Berera, J. C. Collins, Nucl. Phys. B 474 (1996) 183. 

[6] V.A. Petrov, A. V. Prokudin, Eur.Phys.J.C23 (2002)135. 

[7] CDF Collaboration (K. Borras for the collaboration). FERMILAB-CONF-00-141-E, 
Jun 2000; K. Goulianos, talk given in the Xth Blois workshop, 2003, Helsinki, Finland; 
M. Gallinaro, hep-ph/0311192. 

[8] V. A. Khoze, A. D. Martin and M. G. Ryskin, Eur. Phys. J. C 14 (2000) 525; 
Eur. Phys. J. C 21 (2001) 99; V. A. Khoze, hep-ph/0105224; 

[9] V. A. Khoze, A. D. Martin, M. G. Ryskin, Eur. Phys. J. C 26 (2002) 229. 
[10] J.-R. Cudell and O. F. Hernandez, Nucl. Phys. B 471 (1996) 471. 

[II] V. A. Khoze, A. D. Martin, M. G. Ryskin, Eur. Phys. J. C 19 (2001) 477. 



8 



[12] A. De Roeck,V. A. Khoze, A. D. Martin, R. Orava, M. G. Ryskin, Eur. Phys. J. C 
25 (2002) 391. 

[13] V. A. Khoze, A. D. Martin, M. G. Ryskin, Eur. Phys. J. C 23 (2002) 311. 

[14] A. D. Martin, M. G. Ryskin and V. A. Khoze, Phys. Rev. D 56 (1997) 5867. 

[15] A. Bialas and V. Szeremeta, Phys. Lett. B 296 (1992) 191; 
A. Bialas and R. Janik, Z. Phys. C 62 (1994) 487. 

[16] J. Pumplin, Phys. Rev. D 52 (1995) 1477. 

[17] M.G. Albrow and A. Rostovtsev, FERMILAB-PUB-00-173. Aug. 2000. 



9 



Figure captions 

Fig. 1: The full unitarization of the process p + p — > p + X + p. 

Fig. 2: The process p + p-^p + H + p. Absorbtion in the initial and final pp-channels 
is not shown. 

Fig. 3: t-distribution da / dt / a tot of the process p + p — > p + H + p for Higgs boson masses 
100 and 500 GeV. 

Fig. 4: ^-distribution da/d^/a tot of the process p + p — > p + H + p iov M H = 100 GeV. 

Fig. 5: The total cross-section (in fb) of the process p + p — > p + H + p versus Higgs 
boson mass for c gp = 3.5 without Sudakov-like suppression. 

Fig. 6: The total cross-section (in fb) of the process p + p — > p + H + p versus Higgs 
boson mass for c gp = 3.5 with Sudakov-like suppression. 

Fig. 7: The total cross-section (in fb) of the process p + p — > p + H + p versus Higgs 
boson mass for c gp = 2.3 without Sudakov-like suppression. 

Fig. 8: The total cross-section (in fb) of the process p + p — > p + H + p versus Higgs 
boson mass for c gp = 3.3 with Sudakov-like suppression. 
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